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C^ ' Abstract 

In a conformal invariant one-dimensional stochastic model, a certain non-local per- 
^5 I turbation takes the system to a new massless phase of a special kind. The ground-state 

Q ' of the system is an adsorptive state. Part of the finite-size scaling spectrum of the evo- 

CJ . lution Hamiltonian stays unchanged but some levels go exponentially to zero for large 

lattice sizes becoming degenerate with the ground-state. As a consequence one ob- 
serves the appearance of quasistationary states which have a relaxation time which 
grows exponentially with the size of the system. Several initial conditions have singled 
out a quasistationary state which has in the finite-size scaling limit the same properties 
^^ I as the stationary state of the conformal invariant model. 

1 Introduction 

We have recently [1] presented the peak adjusted raise and peel model (PARPM). This is 
'k>( I a one-parameter (denoted by p) extension of the well studied, raise and peel (RPM) one- 

Vh ' dimensional growth model [21 |3]. The latter model is recovered if one takes p = 1. The 

PARPM is defined in the configuration space of Dyck (RSOS) paths on a lattice with L + 1 
sites [L even). The RSOS configurations can be seen as an interface separating a fluid of tiles 
covering a substrate and a rarefied gas of tiles hitting the interface. li h{i) {i = 0,1, ..., L) is 
the height at the site i of an RSOS path, for the substrate one has h{2k) = 0, h{2k + 1) = 1. 
The interface is composed of clusters which touch each other. Depending on the position 
of the hit, the tile can be locally adsorbed (increasing the size of a cluster or fusing two 
clusters) or can trigger a nonlocal desorption, peeling part of a layer of tiles from the surface 
of a cluster. A tile hits a peak (local maximum) with a p dependent probability, and is 
reflected. The other sites are hit with equal probabilities. The effective rates for adsorption 
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and desorption become dependent on the total number of peaks of the configuration, on the 
size L of the system, and on p. 

If the parameter p = 1, the rates are all equal to 1 and independent on the number of 
peaks and the size L of the system. The situation is very different if p 7^ 1. The dependence 
of the rates on the global properties of the configuration can be seen as a process with 
long-range interactions. The larger the value of the parameter p (p > 1), the stronger the 
"long-range" effect is. Configurations with many peaks become more stable. The slowing 
down of configurations with many peaks will lead us to new physics. 

It was shown in [1] that for < p < Pmax, where Pmax = 2(L — 1)/L, in the finite- 
size scaling limit, the properties of the system are p independent. The system is conformal 
invariant (this is the merit of the model) and the stationary states have many well under- 
stood properties. The p dependence of the model appears only in the non-universal sound 
velocity Vs{p) which fixes the time scale. The spectra of the Hamiltonians describing the 
time evolution of the system is given by a known representation of the Virasoro algebra [1]. 
Moreover it was observed that ii p = Pmax the stationary state become an absorbing state, i. 
e., with probability 1 one finds only one configuration. This is a new phase. The absorbing 
state corresponds to the configuration given by the substrate (maximum number of peaks). 
Conformal invariance should be lost. It turns out that the picture is much more complex, 
the PARPM at p = Pmax has fascinating properties. One observes quasistationary states, 
and conformal invariance is broken in an uncanny way. One should keep in mind that one 
deals with a nonlocal model. For p > Pmax some rates become negative and the stochastic 
process is ill defined. 

The present paper deals only with the new phase of the PARPM and is a natural con- 
tinuation of our previous work [1] . The presentation of the model in Section 2 in this paper 
is a mere repetition of Section 2 of [1] . 

Since except for p = 1 the PARPM is not integrable, we have studied its properties 
using Monte Carlo simulations on large systems (up to L = 70000). For the study of the 
spectrum of the evolution Hamiltonian we have done numerical diagonalizations of lattices 
up to L = 18 sites and up to L = 30 for one special case. 

In Section 3, using Monte Carlo simulations we study the time evolution of the system. 
We show that, surprisingly, for moderate system sizes and various initial conditions, after a 
short transient time the system stays practically unchanged for a long relaxation time in a 
quasistationary state (QSS). 

Quasistationary states are observed in systems with long-range or mean-field interactions 
in statistical mechanics and Hamiltonian dynamics. There is a long list of papers on this 
subject and we refer the reader to some reviews O |6]. Typically the time the system spends 
in a QQS increases with the length of the system following a power law (this is not going to 
be seen in our case). We are aware of only one other stochastic model defined on a lattice (the 
ABC model [7]) in which QSS's are seen. In this model in the stationary state translational 
invariance is broken. The sites are occupied in alphabetical order by three blocks of A, B or 
C particles. In the quasistationary states there are more blocks. As we are going to see our 
model is very different. 

In order to understand the origin of the QSS, in Section 4 we do a finite-size scaling 
study of the spectra of the Hamiltonian which gives the time evolution of the system. In 



the unperturbed {p < Pmax) case, the finite-size scahng spectrum is given by the conformal 
dimensions A. They are equal to all non- negative integers numbers except 1 (there is no 
current). The degeneracies are also known [1]. 

It turns out that in the new phase, all the properly scaled energy levels stay unchanged, 
except for one level for each even conformal dimension (this includes the energy-momentum 
tensor which plays a crucial role in conformal invariance and has A = 2 !). We are thus left 
with 

A = 0,3,4,..., (1.1) 

with a degeneracy smaller by one unit for all even values of A starting with A = 4. This is 
not a rigorously proven result but a conjecture based on numerics up to A = 7. 

If we increase the lattice size of the system, all the levels which are not anymore in 
the conformal invariant towers go exponentially to zero. This makes the value A = to 
be infinitely degenerate in the infinite volume limit. One should keep in mind that the 
configuration corresponding to the substrate, which is an absorbing state, has also A = 0. 
For finite volumes, the missing levels are the origin of the QSS. 

In Section 5, in order to make the connection between eigenfunctions and the probability 
distribution functions (PDF) seen in QSS, we mention some properties of intensity matrices 
(the Hamiltonian is one of them) when one of the states is an absorbing state. 

We next derive some properties of the QSS related to the eigenfunction of the energy 
level originally at A = 2 for p < Pmax, and which decreases exponentially to the value A = 
at p = Pmax- This correspondence is possible due to a unique property of the eigenfunction 
of the first exited level of an intensity matrix in the presence of an absorbing state. 

Using Monte Carlo simulations and initial conditions in which the probability to have 
the substrate is taken zero, in Section 6, we have studied the density of contact points and 
the average height in the finite-size scaling limit. If the system is conformal invariant, both 
these quantities are given by precise analytic expressions. It turns out that in the QSS 
state, with very high accuracy, the same expressions describe the data. This result is more 
than surprising. As we have discussed, in the new phase the finite size scaling spectrum of 
the Hamiltonian (which gives the time-like correlation functions) is not the same as in the 
conformal invariant region and one would expect that the space-like correlation functions to 
change too. 

The open questions and our conclusions are presented in Section 7. 

2 Description of the peak adjusted raise and peel 
model 

We consider an open one-dimensional system with L + 1 sites (L even). A Dyck path is a 
special restricted solid-on-solid (RSOS) configuration defined as follows. We attach to each 
site i non-negative integer heights hi which obey RSOS rules: 

hi+i-hi = ±l, ho = hL = (^ = 0,1,...,L-1). (2.1) 



Figure 1: An example of a Dyck path for L = 14. There are four contact points and three 
clusters. The substrate profile is shown in blue. 




Figure 2: Example of a configuration with 4 peaks of the PARPM for L = 18. Depending 
on the position where the tilted tiles reach the interface, several distinct processes occur (see 
the text). 



There are 



Z(L) = L!/(L/2)!(L/2 + l) 
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configurations of this kind. If hj = at the site j one has a contact point. Between two 
consecutive contact points one has a cluster. There are four contact points and three clusters 
in Fig. [U 

A Dyck path is seen as an interface separating a film of tilted tiles deposited on a substrate 
from a rarefied gas of tiles (see Fig. |5]). The stochastic processes in discrete time has two 
steps: 



a) Sequential updating. With a probability P{i) a tile hits the site i = 1, . . . , L — 1 

(Ei ^('0 = !)• In the RPM, P{i) is chosen uniform: P(z) = P = 1/(L - 1). In the PARPM, 
this is not longer the case. For a given configuration c (there are Z{L) of them) with n^ 
peaks. All the peaks are hit with the same probability Rp = p/{L — 1) (p is a non- negative 
parameter), all the other L — 1 — Uc sites are hit with the same probability Qc = Qc/iL — 1)- 
Since 



HcRp + (L - 1 - nc)Qc = 1 
Qc depends on the configuration c and on the parameter p, and we have that 



qc = {L-l -pnc)/{L 



nc), 



1,2, 



ZiL). 



(2.3) 



(2.4) 



b) Effects of a hit by a tile. The consequence of the hit on a configuration is the same 



as in the RPM at the conformal invariant point. Depending of the slope Sj = (/ij+i 

at the site i, the following processes can occur: 

1) Si = and hi > hi^i (tile b in Fig [2]). The tile hits a peak and is reflected. 



■/i.-i)/2 



2) Sj = and h^ < hi_i (tile c in Fig. [2]) . The tile hits a local minimum and is adsorbed 

3) Sj = 1 (tile a in fig. |2]). The tile is refiected after triggering the desorption (hj — > hj — 2) 
of a layer of 6 — 1 tiles from the segment {j = i + 1, . . . , i + b — 1} where hj > hi = hi+b. 

4) Si = —1 (tile d in Fig. |2]). The tile is refiected after triggering the desorption {hj -^ hj — 2) 
of a layer of 6 — 1 tiles belonging to the segment {j = i — b+1, . . . ,i — l} where hj > hi = hi^b. 

The continuous time evolution of a system composed by the states a = 1,2,..., Z{L) 
with probabilities Pa{t) is given by a master equation that can be interpreted as an imaginary 
time Schrodinger equation: 

±P,(t) = -Y^Ha,bm). (2.5) 

The Hamiltonian H is an Z{L) x Z{L) intensity matrix: Ha^b {o- 7^ b) is non positive and 
J2aHa,b = 0. —Ha,h (o 7^ b) is the rate for the transition \b) -^ \a). The ground-state 
wavefunction of the system |0), H\{]) = 0, gives the probabilities in the stationary state: 

\0) = Y.Pa\a), Pa = }imPait). (2.6) 

t— >oo 
a 

In order to go from the discrete time description of the stochastic model to the continuous 
time limit, we take At = 1/(L — 1) and 

Hac = -racqc (c^a), (2.7) 



where rac are the rates of the RPM and qc is given by Eq. (12. 4p . The probabilities Rp 
don't enter in (12. 5 p since in the RPM when a tile hits a peak, the tile is refiected and 
the configuration stays unchanged. Notice that through the QcS the matrix elements of 
the Hamiltonian depend on the size of the system and the numbers of peaks ric of the 
configurations. 

As can be seen from (2.3) and (2.7), for p < 1 the adsorption and desorption are faster 
than at p = 1 and slower for p > 1. The slowing down is extreme for the substrate where 
He = L/2. In this case for the value p = Pmax = 2(L — 1)/L we have Qc = and the substrate 
becomes an absorbing state. 

In a previous paper it was shown that the PARPM is conformal invariant in the domain 
< p < Pmax- The following exact results, which are independent on p, are known for this 
domain. 

The average height for large values of the size L of the system is equal to [51 [5] 

h(L) = - / ' — ln(- sin x)dx + /3^ 0.1056 In L + /3', (2.8) 

n Jf 271 TT 

where /3 and (3' are non universal numbers. 

The density of contact points g{x,L) {x is the distance to the origin), in the finite-size 
scaling limit (x >> 1, L » 1, but x/L fixed) is given by [TU]: 

g{x,L) = Ci-sm{7ix/L)] , (2.9) 
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Figure 3: Special initial conditions for L = 12 (see text). 
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V3 
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r(-l/6) = 0.753149.. 



(2.10) 



The average density of minima and maxima (sites where adsorption doesn't take place) t(L) 
has the asymptotic value: 



lim t{L) = 3/4, 



L— >oo 
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with non-universal corrections (depending on the value of p) of order 1/L. We will use these 
results in Sections 3 and 6. 



3 Quasistationary states at p^^^ 

When we understood that at Pmax one has an absorbing state, and therefore a phase transi- 
tion, we got interested to see how conformal invariance is broken. We expected the system 
to get massive as is the usual case when conformal invariance is broken. It turns out that 
the new phase is a fascinating object. 

We have studied several initial conditions specified by the local heights hi (see Fig. [3]). 
One is the " pyramid^ (PYR): the heights are: 0, 1, 2, 3..., L/2, L/2 — 1, ..., 1, 0. Another is the 
"one denf (OD): the heights being 0, 1, 2, 1, 0, 1, 0, 1, 0, ..., 1, 0. The "two denf (TD) one: 
the heights are: 0, 1, 2, 1, 0, 1, 0, ..., 0, 1, 2, 1, 0. The "all dents" (AD) is defined by the local 
heights: 0, 1, 2, 1, 0, 1, 2, 1, 0, ..., 0, 1, 2, 1, 0. We have chosen these four configurations because 
they are extreme cases. The PYR configuration has only one peak. The OD asymmetric 
configuration has only one peak less than the substrate which has the maximum number 
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Figure 4: p = Pmax- The average height from which one has subtracted the average height of 
the substrate (|), as a function of time for two initial conditions PYR and TD. The lattice 
size L = 96. The averages are obtained taking 6 x 10^ samples. 



of peaks. The symmetric TD configuration has two peaks less than substrate, the AD 
configuration is an intermediate one. 

We have looked, in Monte Carlo simulations, at the average height from which we have 
subtracted the average height of the substrate (equal to 1/2), as a function of time taking 
L = 96, and starting with the PYR and TD initial conditions. The results of our simulations 
are shown in Fig. HI One sees that for each of the two initial conditions one obtains time- 
independent results which are not zero, as one could expect to find in the absorbing state. 
This observation suggests the existence of quasistationary states [5l E]. 

We have also looked at the average density of clusters nci{t)/L from which we have 
subtracted | which is the corresponding quantity for the substrate. The results for the two 
initial conditions PYR and TD are shown in Fig. |5l One obtains similar results as those 
shown in Figure HI 

Comparing the average heights and the density of clusters for the two initial conditions, 
suggest that one has more clusters (therefore lower values of the average height) for the TD 
initial condition compared to the PYR initial condition. Therefore for L = 96 the two QSS 
are different. 

Let us make an observation. If we assume that in the QSS one has the same density of 
clusters as in the stationary distribution observed in the < p < Pmax case [2], for which 
one has an analytical expression {rid/L = T{l/3)\/3/27rL^^^), one obtains the value —0.339 
for the quantity shown in Fig. [5l This value is closed to the value seen for the PYR initial 
condition. This observation will play in important role in understanding the QSS. 

In order to see how the QSS appeared, we looked for another quantity which was ex- 
tensively studied in the PARPM [T]. This is the average density of sites where one has a 
maximum or a minimum in a given configuration T^i^x{L,t) for p = Pmax- In the PARPM 
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Figure 5: The density of clusters - 1/2 (the density of clusters in the substrate) as a function 
of time. Same conditions as in Figure 3. 



and < p < Pmax! in the large L limit one has r(L, oo) = r(L) = 0.75 (for the substrate 
TiL) = 1). 

If one starts with the TD configuration and looks at time variation of the quantity 
1 ~ Tma.xiL,t) for different system sizes one obtains the results shown in Fig. El We can 
see that for small values of L = 46 — 60 one has an exponential fall-off with time. Then 
an almost linear decrease in time for the values L = 62 — 76. For L = 92 and 96 one sees 
practically no time variation and a value rmax(-^, t) ~ 0.77 very close to the value 0.75 seen in 
the stationary state in the conformal invariant domain of the model. This makes us suspect 
the following behavior of T^ax{L,t) in the QSS: 



rmax(^,t) = A(L)exp(-E(L)t), 



(3.1) 



where E{L) decreases exponentially with L and A{L) increases smoothly with L to the value 
0.25 observed for the stationary states in the < p < Pmax domain. In Fig. [TJ we show for 
L = 50 the data undistinguished from a fit: 

1 - rmax(50,t) = 0.217exp(-0.0000522t). (3.2) 

Notice the very small value of E and the fact that y4(50) is not far away from the value 0.25. 

We sum up our observations; although for p = Pma.^ we expected the system to relax in the 
absorbing state, we observed the existence of states with very long relaxation times. For the 
lattice sizes presented above, the QSS depend slightly on the initial conditions. Surprisingly, 
the density of clusters and of maxima and minima in the QSS have values closed to those 
observed in the stationary states for < p < Pmax- Actually for the PYR initial condition 
the results coincide. In the next three Sections we will explain these observations. 
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Figure 6: p = Pmax- Average density of minima and maxima T^i^^{L,t) subtracted from 
1, as a function of time for different lattices L. The initial condition OD was chosen. The 
averages are obtained by taking 6 x 10^ samples. 
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Figure 7: The average density of minima and maxima t(L) as a function of time for L = 50 
together with the fit (3.2). The initial configuration is OD and the averages are obtained by 
taking 6 x 10^ samples. 



4 The spectrum and wavefunctions of the Hamilto- 
nian. How QSS occur at p^ 



max* 



In order to understand the origin of the QSS we have studied the spectrum of the Hamiltonian 
which gives the time evolution of the system (see Eq. (2.5)). For < p < Pmax, where we 
have conformal invariance, we have checked [T] that in the finite-size scahng hmit: 

hm E,{L) = 7rt;,A,/L, z = 0, 1, 2, . . . , (4.1) 

where Eq = 0, Aj are the scahng dimensions, and the sound velocity has the expression: 

Vs{p) = {l-3{p-l)/5)3V3/2. (4.2) 

Notice that the velocity decreases when p > 1 since, as described in Section 2, the transition 
rates are smaller. The scaling dimensions are given by the partition function [1]: 



oo 



Z{q) = T.<l^^ = {l-q)m^~^'r'- (4.3) 

i=0 n=l 

We give the first values of A^'s together with the corresponding degeneracies (rf^'s): 

A = 0(1), 2(1), 3(1), 4(2), 5(2), 6(4), 7(4), .... (4.4) 



^inax- 



We will check if these values will be seen also for p = Pn 

In order to estimate the values of the Aj's, we have taken L = 18 (this is not a small 
lattice!) and diagonalized numerically the Hamiltonian for various values of p. The results 
are shown in Fig. [8] where the first 11 levels are seen (the ground-state energy Eq is equal to 
zero). The remaining 10 levels should correspond roughly (see (4.4)) to A = 2,3,4,5 and 6. 
We can see that for p = 1 where the model is integrable, this is indeed the case. The levels 
cluster in the right places. When p increases, one notices that the properly scaled Ei after 
a smooth behavior up to p ~ 0.9, decreases rapidly for p = Pmax (we have used for fs(pmax) 
the value given by (4.2) for p = 2). Using Monte Carlo simulations, we have checked [1] that 
for p < pmax the small decrease with p of Ei is a finite-size effect and therefore what one sees 
in the figure is a crossover effect. One can also see that for increased values of p, E^ crosses 
£^3 and that Eg crosses all the levels Eg, . . . jE^. Except for the three levels Ei, E4 and Eg 
which decrease dramatically for p = Pmax, the other levels have the same finite-size behavior 
as those in the conformal invariant domain [p < Pmax)- This suggest that the three levels 
mentioned above might be related to QSS. We now proceed to a detailed analysis of these 
observations. 

Using different lattice sizes we have computed Ei as a function of L up to L = 30. For 
this calculation we could study larger lattices due to a special property of the Hamiltonian 
aX p = Praia- As wc are going to show in section 5, the eigenlevel corresponding to Ei is 
the ground state energy of a reduced matrix defined in a basis where the absorbing state is 
absent. In this case, by using the power method we were able to calculate Ei up to L = 30. 
The results can be seen in Fig. O Using the two points corresponding to the largest lattice 
sizes, one obtains: 

Ei{L) = 0.912 exp(-0.206L). (4.5) 

10 




(p-1)/(p^..-1) 



Figure 8: Estimates A^^^ = LEn/nVs^p) of the scaling dimensions A„ for different values of 
p, and for the lattice size L = 18. The estimates corresponding to the first 11 energy levels 
are shown. The values of Vs{p) were obtained by (4.2). 
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Figure 9: ln(£'i) as a function of the lattice size L ioi p = Pmax- The red line is a guide for 
the eyes, and is obtained from a fit where the lattice sizes L = 16 — 30 were used. 
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Figure 10: Estimates of —d/ln{L){\n{Ei{L))} as a function of 1/L for different lattice sizes 

(P = Pmax)- 



To be sure that a power law behavior {E ^ L~™) is excluded, we have estimated the 
derivative — (i/(iln(L){ln(i^^i(L))}. This quantity should reach a constant for large values 
of L if Ei{L) is power behaved but should diverge linearly in the case of an exponential 
behavior like (4.5). As seen from Fig. [TOl the exponential behavior (4.5) is correct. 

A similar analysis of E^{L) but up to L = 18 only, gives a similar result: 



E^{L) = 2.41exp(-0.10L). 



(4.6) 



We have not looked at Eg{L) but we expect again an exponential fall-off. We conclude that 
three energy levels have an exponential fall-off and that they can be related to QSS. This 
is going to be shown to be the case in Section 5. We proceed by looking at the remaining 
levels. 

The data suggest that at p = Pmax, the energy levels that do not go exponentially to zero 
have the same finite-size scaling behavior as those in the conformal invariant domain. This 
would imply that instead of (4.4) we would have 



A = 0(l),3(l),4(l),5(2),6(3),... 



(4.7) 



If confirmed, this would lead us to a strange picture since the scaling dimension A = 2 does 
not appear. This dimension corresponds to the energy-momentum, and therefore conformal 
invariance couldn't apply and we could not explain the finite-size behavior of the remaining 
levels. What can go wrong in our picture? One possibility is that the finite-size scaling of 
the levels doesn't satisfy (4.1). 

We have computed E2{L) up to L = 18. A fit to the data gives A3 = 3.05 in agreement 
with what should be expected. Similarly examining E^^L) one finds A4 = 3.95 also as 
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Figure 11: Ratio among the two lowest eigenenergies, at p = Pmax,that does not vanish 
exponentially. The data are plotted as a function of 1/L, for lattice sizes L = 8 — 18. A 
cubic fit was done by using the five largest lattice sizes. 

expected. These estimates were found assuming that fs(pmax) is given by Eq. (4.2). Can 
we get A = 2 changing the sound velocity such that E^{L) gives A = 2, E^{L) gives A = 3, 
Eq{L) gives A = 4, . . .? We have computed the ratios E^{L) / E^{L) and Eq{L) / E^{L) as a 
function of L. One should obtain 4/3 respectively 5/3 if one had (4.3) and 3/2 respectively 
2 if the energy momentum tensor would be present. In Figs. [11] and [12] we show these ratios 
as functions of 1/L. Cubic fits give the values 1.32 respectively 1.64. We conclude that (4.7) 
is most probably correct. We have also checked that there are no energy crossings for the 
levels which cluster around A = 7. 

The analysis of the energy levels (compare (4.4) with (4.7)) suggests that aX p = Pmax 
the partition function (4.3) changes in the following way: the degeneracy at each even value 
of A decreases with one unit. Each energy level which left the Virasoro representation at 
non-zero even values of A moves to A = which becomes infinitely degenerate. This opens a 
problem in the representation theory of the algebra which might be solvable since the central 
charge is c = 0. In Section 6 we are going to learn more about this puzzle. 

5 From eigenfunctions of the Hamiltonian to QSS 

We have seen in the last section that some eigenvalues of the Hamiltonian vanish exponen- 
tially. Here we will show what is their connection to QSS. In order to do so, we first prove 
a special property of Hamiltonians in the presence of an absorbing state. 

We denote the vector space corresponding to n+1 configurations by |0), |i) (i = 1, 2, ..., n), 
in which we have chosen |0) to be the absorbing state. The Hamiltonian has the following 
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Figure 12: Ratio among the third and second lowest eigenenergies, at p = p^ax^that does not 
vanish exponentially. The data are plotted as a function of 1/L, for lattice sizes L = 10 — 18. 
A cubic fit was done by using the four largest lattice sizes. 



properties: 



Hifl = 0, Hofi = 0, Hij < 0, 
Ho,, + J2H^,, = (j = l,2,...n). 



(5.1) 
(5.2) 



Assume that E^. {k = l,...,n) is a non- vanishing eigenvalue of H with an eigenvector 
{Uo )l/i \---^yn'')- We can show that the sum of the components of any eigenvector is 
equal to zero: 

2/f + E2/f = (fc = l,2,..,n). (5.3) 



Using (15. ip and (15. 2p we have: 



,(fc) 



E,yr = E Ho.yr = - E H,,iyr = -Eu E v^ 



(fc) 



,(fc) 

3 ' 



(5.4) 



j.« 



from which the identity (5.3) follows. 

We consider now the reduced matrix H[- = Hij {i,j = 1, 2, ..., n) (the configuration |0) is 
taken out). Let Ei be the lowest non- vanishing eigenvalue, from (15.10 and using the Perron 
Frobenius theorem, we get: 

yf ^ > 0, (5.5) 

and using (5.3), yg < 0. From the same theorem we also learn that Ei is the unique 
eigenvalue for which (15. 5p occurs. For the other eigenvalues, at least one component of the 
wavefunction is negative, i. e. Eq. (5.5) is not valid for yl (k > 1). 
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The solutions of the differential equations f l2.5p are 

Po{t) = I + Y,A,yi!'^ exp{-Ekt), 

k 

Pi{t) = J2 ^fcl/f exp(-Efct). (5.6) 

k 

The n constants A^. are determined from the initial conditions. 

At large values of L and t, the exponentially falling energies give the major contributions 
to Po(t) and Pi(t). Among them, Ei plays a special role. Note only is Ei the smallest 
energy, but the components of its eigenfunction are positive (5.5). This implies that for a 
large range of L and t (both of them large), one can keep only the term with fc = 1 in the 
sums appearing in (5.6). The situation is different at very large values of L when all the 
exponentially falling energies are practically equal to zero and more terms can appear in 
(5.6). Independent of the initial conditions, the term with k = 1 has to be present in the 
sums (5.6) in order to assure the positivity of the probabilities Po(t) and Pi(t) since for the 
other eigenf unctions (5.5) is not valid. For example, the eigenfunctions corresponding to the 
levels Ei{L) and Eg{L) of Fig. [8] have components with both signs in the sum giving Pi{t). 

We are now in the position to show how exponentially falling energies like Ei[L) (see 
(14. 5p ) can be the origin of QSS. For our discussion we will assume that the term with k = 1 



alone is present in (5.6). It turns out that this assumption will help to explain all the results 
obtained for the Monte Carlo simulations. The probability vector |P(t,L)) becomes: 

|P(t,L)) = [l-a(L)exp(-Ei(L)t)]|0) + a(L)5]p,(L)|z)exp(-Ei(L)t), (5.7) 

i 

where we have used (15. 3p and (15.50 

a{L) = A,Y.yf\L); ML)=yi'\L)/j:yl'\L). (5.8) 

i i 

Note that Pi{L) gives the probability to find the system in configuration \i) if the system is 
in the stationary state of the reduced Hamiltonian H', that acts in the vector space where 
the substrate is absent. Thus, Equation (15. 7p explains the occurrence of QSS. If L is large 
enough, Ei{L) is negligible, one finds: 

\P{L)) = (1 - a(L))|0) + a{L)\PUL)). (5.9) 

Here a{L) depends on the initial conditions and it is not equal to zero. |Pns(-^)) is a proba- 
bility distribution function of configurations in which the substrate is not present. Equation 
( 15. 9 p describes therefore the QSS. Using Monte Carlo simulations we have measured the 
probability to find the system in the substrate. For the OD initial condition (see Section 3 
for the definition), one finds: 

1 - a(L) ~ 4.5/L. (5.10) 

For the PYR initial condition one finds 1 —a{L) = 10~^ for L = 96 (we didn't compute a{L) 
for the TD and AD initial conditions). This implies that for large values of L one has 

|P(L)) = |P,,(L)), (5.11) 

and the substrate does not occur in the QSS. We postpone the discussion of the results 
shown in Fig. [6] up to the next Section. 
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Figure 13: The average height h{L) as a function of ln(L) in the QSS. L = 1024, 2048, 4096 
and 8192. We have used PYR as an initial distribution. 

6 The quasistationary states 

In this section we are going to identify the quasistationary states observed in the Monte Carlo 
simulations (see Figs. HUT]) and find a puzzle. The basis of our identification are Eqs. (5.7), 
(5.9) and (5.11). 

We have taken large lattices and started our simulations with the PYR initial condition. 
We first looked at the average height h{L). The results are shown in Fig. [T31 The data can 
be fitted by a straight line. Taking only the two largest lattice sizes we obtain: 



h{L) = 0.1068 ln(L) + 0.042. 



(6.1) 



In astonishing agreement with the expression (2.8) derived [HI [9] assuming conformal in- 
variance and used to describe the data for < p < Pmax- Since for the large lattices we 
have considered, we can use (5.11) and this would imply that |Pns(-^)) is given by the same 
function as the one seen away from Pmax- 

If confirmed, this result would be surprising because as we have discussed in Section 4, 
as compared to the conformal invariant domain, at p = Pmax the finite-size spectrum of the 
Hamiltonian is a "mutilated" one. It lacks not only the energy momentum tensor but other 
levels too. On the other hand that space-like correlation functions look to be unaffected. 

In Fig. [13] we show for the AD initial condition, the density of contact points g{x, L) 
for various lattice sizes divided by the finite-size scaling distribution (2.9) in the QSS. The 
coincidence of the Pmax data and the expectation coming from conformal invariance in the 
QSS is astonishing (for the large lattices considered here a{L) is practically equal to 1). 

From now on, we will assume that for large lattices the QSS the correlators are those of 
the conformal invariant model (RPM) and will try to explain the results described in Section 
3. 



16 



15 



>< 

>f^ 



1.1 



1.05 




■| 0.95 



o» 



0.9 



0.85 



L=512 
L=1024 
L=2048 
L=4096 



0.2 



0.4 0.6 

x/L 



0.8 



Figure 14: The density of contact points g{x, L) as a function of x/L at Pmax? divided by 
(2.9). L = 1024,2048,4096 and 8192. Initial condition AD. The average are obtained from 
3500 samples running 10^ mcsteps. 



We use Eq. (5.9) and the observation that for the initial condition PYR a{L) = 1, taking 
a(96) = 0.975 for the TD initial condition we obtain the results for TD presented in both 
Figs, m and [5] (we don't have at hand an independent measurement of a{L) for the TD initial 
condition) . 

We proceed by looking for a description of the data presented in Fig. El where the time 
dependence of the quantity 1 — Tma.x{L,t), taking OD initial condition, is shown. Traax{L,t) 
is the average density of minima and maxima which is equal to 1 for the substrate. We use 
Eqs. (5.7) and (5.10) to get: 



1 - r„,ax(i^, t) = [-4.5/L + 1 - (1 - 4.5/L)r(L)] exp(-Ei(L)t). 



(6.2) 



The function r(L) is known analytically only for p = 1 [2]. It is equal to 3/4 in the large 
L limit and has non- universal corrections of order 1/L which are unknown for p = Pmax- In 
principle these corrections can be determined by looking at the QSS obtained using PYR 
initial conditions. We can however use the result obtained in [1] for p = 1.99 where we found 
t{L) ^ 0.75 — 0.3/L and use it in (6.2) to get a fair description of the data. We obtain 

l-r.,,.(L,t) ^ (0.25-0.8/L)exp(-^i(L)t). (6.3) 

The value of -E'i(L) can be estimated from (4.5). In Fig.fT^we use this function, together with 
the predicted values of Ei{L) given in (4.5), to compute the time dependence of 1 — rjnax(-^, t) 
for the same lattice sizes used in Fig. [61 We can see that the overall behavior of Fig. [151 and 
Fig. [6l generated by the Monte carlo simulations, are the same. 



17 



0.25 


- 














L-96 


-1 




P 


s^ 


^ 


^= 


=z 


_ 





L^76 


- 


0.2 


1 


\\ 
\ \ 

w 


\ 


^...^^ 


^ 


^ 




^^r^ 


--^ 


J 0.15 


1\ 








\ 


\ 




-^^^^^^~ 


•-..^^ 


H 


1 


\ 


\ 








^ 


"t««r"~~^ 


-^- 


" 0.1 


:\ 


^\ 


\ 






^ 


^ 
^ 


v^\ 


- 


0.05 


- 


\ 


\ 


\ 


L^ 


^ 




"^^-...^^J"'^^^^ 


--- 




- 


l\ 


^ 


L^ 


-^ 






— h- III, 


--- 



20000 40000 



60000 
t 



80000 



le+05 1.2e+05 



Figure 15: The average density of minima and maxima T^i^^{L,t){L,t) subtracted from 1, 
as a function of time for the same lattice sizes L, as in Fig El The initial condition is OD. 
The data were generated from the prediction (6.3), using Ei given by (4.5) 

7 Conclusions 



The parameter p which enters in the definition of peak adjusted raise and peel stochastic 
model determines a domain < p < Pmax in which one has conformal invariance. In the finite- 
size scaling regime all the properties are independent on p which appears only in the sound 
velocity Vs{p) which fixes the time scale. If p > 1, global properties of the configurations 
(the number of peaks) and the size of the system L, enter in the rates. The larger is the 
number of peaks of a configuration, the smaller are the rates to escape the configuration. 
As a result, at p = Pmax = 2(L — 1)/L the configuration with the largest number of peaks 
which is the substrate, becomes an absorbtive state for any size of the system. Since there 
are no fluctuations in the ground-state of the system, we expect conformal invariance to be 
lost and get into a massive phase. This is not the case and a fascinating phenomenon takes 
place. 

The new features of the model at p = Pmax are encoded in Figures H] and [Hi If the 
evolution of the system starts with a given initial configuration, instead of moving fast to 
the absorbing state, the system gets stuck in another configuration and the relaxation time 
grows exponentially with the size of the system. This is a quasistationary state. In Figure 
[8] we look at the spectrum of the Hamiltonian and follow the change of the scaled energies 
with increasing values of the parameter p. From the eleven energy levels shown in the figure, 
eight of them are those of the conformal invariant region {0 < p < Pmax)- Three of them go 
exponentially to zero for large values of L. 

In this paper we have tried to clarify this phenomenon. Based, unfortunately, on numer- 
ics, the following picture emerges. The wavefunction corresponding to the first excited state 
whose energy vanishes exponentially, gives the probability distribution function of the QSS. 
This is due to a peculiar property of the first exited state of any stochastic Hamiltonian in 



the presence of an absorbing state. For finite values of L, the QSS depends on the initial 
conditions but in the thermodynamical limit, it becomes independent of them. Unexpect- 
edly, the finite-size scaling properties of the QSS are identical to the one observed for the 
ground-state in the conformal invariant domain. 

This picture describes a strange way to break conformal invariance. Part of the spectrum 
stays unchanged (not the scaling dimension of the energy-momentum tensor!) and another 
part sinks in the vacuum which becomes probably infinite degenerate. We have not a clue 
how to explain this observation. At the same time, the space-like correlations seen in the 
QSS are those of the conformal domain. This is obviously an unexplained paradox. The 
study of space-time correlation functions should shade light in this problem. We plan to do 
it in the future. 

We have studied only four initial conditions and found essentially only one QSS. This 
QSS could be related to the eigenfunction corresponding to the first exited state. It is most 
probable that we have missed other QSS which should be related to some linear combination 
of eigenfunctions corresponding to the remaining exponentially decaying energies. 

The search for QSS should continue in some extensions of the PARPM. One can look at 
the fate of defects like those studied in [TT] when the rates are adjusted to the number of 
peaks. Another interesting avenue is to study the effect of changing the rates depending on 
the number of peaks in the extension [2], [3] of the raise and peel model in which the rates of 
adsorption and desorption are not equal. 
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